Glassy behavior in systems with Kac-type step-function interaction 
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We study a system with a weak, long-range repulsive Kac-type step-function interaction within 
the framework of a replicated effective <p 4 -theory. The occurrence of extensive configurational en- 
tropy, or an exponentially large number of metastable minima in the free energy (characteristic of 
a glassy state), is demonstrated. The underlying mechanism of mesoscopic patterning and defect 
organizations is discussed. 
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Competing interactions on different length scales cause 
in many cases the emergence of an intermediate length 
scale where new structures and inhomogeneities are 
formed. Examples are stripe formation in doped Mott 
insulators M, bubbles of electronic states of high Lan- 
dau levels in quantum Hall systems 0] , domains in mag- 
netic multilayer compounds H, and mesoscopic struc- 
tures formed in self-assembly systems M . These systems 
typically exhibit a multi-time-scale dynamics similar to 
the relaxation found in glasses. The glassy behavior and 
the diverging relaxation time are believed to be the result 
of the competition between the interactions with different 
characteristic length scales |^,||]; for example the macro- 
scopic phase separation is frustrated by the competing 
long-range interactions |5|,fjJ . Glassiness then arises spon- 
taneously in the absence of extrinsic disorder due to self- 
generated randomness. 

While there are various different scenarios for glassy 
behavior such as the kinetic constraint where the diverg- 
ing relaxation time is purely of dynamic origin and occurs 
in a system with trivial equilibrium properties, the cen- 
tral theme of this work is based upon the random first 
order transition where glassiness is attributed to an 
exponentially large number of metastable states, origi- 
nally emphasized by Kauzmann |^]. The fact that con- 
figurational entropy is needed for slow motions in glasses 
was first described by Gibbs and DiMarzio flQfl . Be- 
low a crossover temperature Ta, an energy landscape- 
dominated, "viscous" long-time relaxation sets in due 
to an exponentially large number of metastable states, 
Af, i.e., the configurational entropy, S c = InA/", becomes 
extensive. This cross-over temperature Ta is often as- 
sociated with the mode-coupling temperature at which 
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the relaxation time or the viscosity exhibits a power 
law divergence \T — Ta| -7 , within mode-coupling the- 
ory [pd| . Activation processes, which are neglected in 
mode-coupling theories, soften the sharp transition into a 
crossover where for T < Ta free energy barriers and thus 
transition rates between the metastable states remain fi- 
nite. The configurational entropy decreases with decreas- 
ing temperature, and becomes negative for T < Tk- A 
continuous (random first order) transition, the "ideal" 
glass transition, occurs at Tk to avoid S c < 0. The 
Kauzmann temperature Tk is the temperature at which 
the genuine thermodynamic glass transition is expected, 
whereas the experimentally observed glass transition oc- 
curs at T g > Tk which depends on the cooling rate. 

In Ref. p], it was shown quite generally within the 
framework of a replicated <^ 4 -theory that the competi- 
tion between short-range forces (favoring phase separa- 
tion) and long-range Coulomb interaction leads to an ex- 
ponentially large number of metastable states and self- 
generated glassiness. The large phase space of fluctua- 
tions, which can lead to a fluctuation induced first or- 
der transition [p^j , was shown to, alternatively, drive the 
system into an amorphous state, the stripe glass. It is 
of particular interest to explore to what extent explicit 
competition or frustration is necessary to cause a glassy 
state and what kind of interactions support such a state. 
A particularly interesting potential is the Kac-type step 
potential 

V(x) = aV<X7z), (1) 

with 4>{y) = 1 for y < 1 and zero otherwise. Here, 7 con- 
trols the amplitude and range of the potential, whereas 
the inverse length a characterizes the integral strength 
J d D xV(x) oc a 2 . In the van der Waals limit, 7^0 
after the thermodynamic limit, it is known that a sys- 
tem of particles interacting through a potential given by 
Eq. (^) can be described exactly by a mean-field theory 
p3[ , where a is related to the long-range force of the van 
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der Waals theory. For the particular case of a step poten- 
tial, mean-field theory predicts a spinodal. This model 
system is of particular interest since it has been used to 
study the glass formation and crystallization processes 
jl4]]. Using Monte Carlo simulations, the appearance of 
many metastable amorphous "clump" configurations was 
demonstrated in Ref. E3|. However, it is difficult, if not 
impossible, to enumerate the metastable states and to de- 
termine the dependence of the number of local minima, 
as a function of the system size, in a Monte Carlo study. 
We find that an analysis of a replicated <p 4 -theory along 
the lines of Ref. [0] is useful as an alternative strategy. 

In this paper we demonstrate that the Kac-type step 
function interaction, Eq. (|l|), indeed causes the emer- 
gence of an exponentially large number of metastable 
states and a self- generated glassy state. In the limit of 
small but finite 7 no actual frustration between differ- 
ent interactions is necessary, which is different from the 
model studied in Ref. H|). We demonstrate that glassi- 
ness is due to multiple configurations of self-arrested de- 
fects, and occurs, similar to Ref. 0, once the correlation 
length of the system is slightly larger than the length 
scale of mesoscopically modulated structures. The con- 
figurational entropy behaves for D = 3 like S c /V oc 7 3 , 
where, V is the volume of the system. A model of glassy 
behavior with the Kac-type long-range interaction, called 
a van der Waals glass, was introduced and extensively de- 
veloped, including a study of its dynamics, in Ref. Jl5| , 
where it is argued that proximity to the mean-field spin- 
odal provides a long correlation length in addition to the 
Kac potential range, and thus leads to frustration and 
nonzero configurational entropy. This will be verified 
here using the replica approach. 

We start from the model Hamiltonian 



H 



+ [ d 3 x { d 3 x'<p(x)<p(x')V(x-x'), 
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with V(x) of Eq. (^). The usual equilibrium free energy 
J- = —T In Z is an outcome of an unconditional average 
over the entire configuration space. It does not permit 
the detection of local minima of the free energy in the 
configuration space. In Ref. |l6f| , a replica approach is 
proposed to overcome this limitation and allows us to 
probe the number of metastable states in models with 
the additional long-range interaction V(x). We adopt 
the self-consistent screening approximation to the repli- 
cated (^-theory developed in Ref. ](|, where the free 
energy J-(m) of the replicated Hamiltonian is given in 
terms of the regular correlation function G(q) and the 
correlation function F(q) = ( ( p a (q)f b (-~q)) between the 
fields in different replicas. Here, a 7^ b are the replica 
indices and m is the number of replicas. The quan- 
tity F(q) equals the Edwards- Anderson order-parameter, 



lhnt^oo{tp(—q,t)ip(q,0)}, which becomes nonzero in the 
glassy state, or when ergodicity is broken. 

Within the self-consistent screening approximation, 
the free energy J-(m) can be expressed as 



J-(m) 



T 



(trine -1 + trlnX>~ 1 ), 
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and the configurational entropy S c is related to ^{m) 



by S c = ^ ^ m ' . Q is the generalized correlation 

m=l 

function matrix, defined by Q = (G — F)I + F~E, with the 
constant matrices I a & = <5 a & and E a f, = 1. The symbol 
tr in Eq. (||) includes the trace of the replica and the 
momentum integral. The momentum dependence of the 
matrices is implied here but omitted for convenience. The 
matrix V is related to Q via the relation V^ 1 = (uT)^ 1 + 
n, where n is the generalized polarization matrix defined 
asn = (G ®G - F ® F)I + (F ® F)E, where ® denotes 
convolution in Fourier space. The replicated Schwinger- 



Dyson equation can be written as Q 1 
£ is the self-energy and 



G (q) = 



1 



r 



v( q y 



G where 



(4) 



with the renormalized mass r — r + j n-ir) 3 ^(?) 
denotes the Fourier transform of the potential V(x). 
Within the self-consistent screening approximation the 
self-energy has diagonal elements T,q — 2G ® Dq and 
off-diagonal elements 'Ep = 2F ® Dp in replica space, 
where Dq and Dp being, respectively the diagonal and 
off-diagonal elements of the matrix T>. These equations 
form a closed set of self-consistent equations which en- 
able us to solve for G and F, and then determine the 
configurational entropy via 
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Here, S[x] = —x — ln(I — x). It is obvious from the above 
equation that F ^ implies S c ^ 0. 

In the limit of small and large 7 it is possible to make 
analytic progress. For large 7, i.e. short-range poten- 
tials the gradient term in Eq. (^) dominates and V(q) 
can be neglected. No glassy state with finite Edwards- 
Anderson parameter results. The situation is more inter- 
esting in the limit of small 7, i.e. for long-range interac- 
tions. Using V(q) = a 2 ip(q/j) with ip(z) — 47r(sin(z) — 
z cos(z))/z 3 , it follows for ]«a that the gradient term 
in Eq. (Q) can be neglected compared to the Kac-type in- 
teraction. The long wavelength behavior is dominated by 
the long-range interaction. The short-range interaction 
becomes effectively local and has no characteristic length 
scale anymore. The correlations are dominated by wave 
vectors which minimize V(q). Since if)(z) is minimal (and 
negative) for z = zq = 5.76, the dominant peak in the 
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FIG. 1. Potential used in the numerical calculations (solid 
line) and the step-function potential with 7 = 0.03 and a = 1 
(dashed line). 
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FIG. 2. Configurational entropy S c (within a replica the- 
ory) as a function of temperature for ro = 0.1494 and 
u = 1.79. 



correlation function occurs at go = z ol, independent of 
a: 



Go (q) 



Z 



e- 2 + (|q|-9 r 
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Here Z = -(7/a) 2 is the weight of the peak with width 
characterized by £~ 2 = Z{r + ca 2 ), where c = ip"(zo) « 
0.361. This expansion around q = qo elucidates the cor- 
respondence with the analysis performed in Ref. ||. A 
calculation along the lines of Ref. || determines T\ as 
the temperature where the ratio of the correlation length 
£ and the modulation length Iq — 2n/qo becomes larger 
than 2. Inserting Eq. (JsJ) into Eq. (||) yields for the con- 
figurational entropy S c — VCj 3 , with C « 6.81 x 10~ 3 . 
The longer range the interaction, the smaller is the num- 
ber of metastable states per unit volume. As expected, 
glassiness disappears in the van der Waals limit 7^0 

0- 

The set of self-consistent equations is solved numeri- 
cally, without assuming Eq. (^J) . A slightly modified Kac- 
potential (see Fig. |l|) has been used for convenience. The 
major conclusions of this paper seem not affected by this 
modification. In Fig. @, we show the temperature depen- 
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FIG. 3. Instantaneous correlation function G(q) (solid line) 
and long-time correlation function F(q) (dashed line) for 
T = 1, ro = 0.1494, and u = 1.79. The arrow points at 
the position of the first secondary peak. 



dence of S c for r — 0.1494 > and u = 1.79 (we use 
units where the upper cutoff of the momentum integra- 
tion is unity) . The behavior of S c matches identically the 
entropy crisis scenario of the random first order transition 
theory |8|. The mode-coupling temperature Ta and the 
Kauzmann temperature Tk can be unambiguously iden- 
tified. Furthermore, we show in Fig. ^ the correlation 
functions G(q) and F(q). In addition to the pronounced 
peak at go we can also see higher order structures in the 
instantaneous correlation function G(q), which become 
strongly suppressed in the long-time correlation function 
F(q). The typical length scale which determines the sup- 
pression of long-time correlations for larger q is the defect 
wandering length A which results from the off-diagonal 
self-energy in replica space Sp(go) = — (j) |§- 

Our analysis allows us to determine what kind of inter- 
actions can cause a glassy state. For small 7 the short- 
range term will be effectively local. Competing inter- 
actions of different length scales can satisfy the glass- 
forming condition. Their presence is, however, not nec- 
essary. We have just demonstrated that nonzero config- 
urational entropy can arise even if there is no explicit 
competition between the interaction length scales. In 
this case, not only does V(q) possess a minimum at go, 
but also V(qo) is negative. The system exhibits a spin- 
odal when ro > 0, i.e. for a purely repulsive interaction, 
in distinction to the model discussed in Ref. (|||. The 
glass-forming criterion £ > 21$ can be fulfilled even in 
the absence of the gradient term. We recall that prox- 
imity to a mean-field spinodal provides a large £ in addi- 
tion to the Kac potential range, and thus self-consistently 
leads to frustration. One can in fact "map" the van der 
Waals glass |TJ[ into the stripe glass It is thus pos- 
sible to predict if glassiness occurs in a system interact- 
ing through a mesoscale potential V(x) . For example, 
in the case of a Kac-type Gaussian potential [Q with 
4>{x) — exp(— x 2 ), the Fourier transform of this potential 
is monotonic. A modulation length ^o at the mesoscale 
does not occur and no glassy behavior is anticipated in 
the absence of the short-range gradient term. This agrees 
with the conclusion reached in Ref. pj . Both the gradi- 
ent term in the Hamiltonian and ro < are necessary in 
order to find similar glassiness at mesoscales. 

From these considerations it seems to follow that the 
free energy landscape of metastable states which causes 
glassiness originates from configurations of mesoscopic 
defects. The origin of the glassiness has been discussed 
in detail in Ref. H in terms of a novel length scale, the 
defect-wandering length A, which depends on £ and Iq. In 
the parameter region where the system acquires modu- 
lation at Iq, the mesoscopically ordered state constitutes 
the global minimum of the free energy J 7 , or the ground 
state. With respect to this ordered array of mesoscopic 
structures, excitations such as 'dislocations' are termed 
defects. The defect-wandering length A is the distance 
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on which a defect can move freely in the lattice of the 
mesoscopic structure. When A < 2Z /3 0], the defects 
are pinned by the underlying lattice of the mesoscopic 
structure. A distribution of such pinned defects becomes 
a local minimum of T . The organizations of defects 
should be responsible for the exponentially large number 
of metastable states M — expS^, where S c is extensive. 
This picture is consistent with the findings of Ref. 
that the FCC arrangement of the clumps has lower free 
energy than those of all frozen amorphous clump config- 
urations. In this case, the FCC clump configuration is 
an ordered array of mesoscopic spherical structures. The 
amorphous clump phases are the organizations of defects 
about the ordered lattice of these mesoscale spherical 
structures. 

Another interesting question is whether glassiness is 
also possible with a microscopic Iq on the order of the 
hard core radius of the atoms. The defect-pinning pic- 
ture for the metastable configurations described above is 
intuitively clear and consistent with the 3D models where 
l is mesoscopic. Nothing in the argument forbids its ap- 
plication to a circumstance where Iq is the microscopic 
lattice constant. To study this in more detail, we have in- 
vestigated a ID 95 4 -model on a lattice in which the frozen 
configurations of kinks are expected to give a large num- 
ber of metastable states. This occurs in the parameter 
regime where the continuous variable tp can be mapped 
to the ID Ising model with nearest-neighbor interaction, 
which has been shown to have S c — . A crude esti- 
mate of the number of low- lying metastable states reflects 
algebraic, instead of exponential size-dependence, leading 
to limy^oo S c /V — 0. While the effect of dimensionality 
is not clear, it is thought not to change the system size 
dependence of the number of low-lying metastable states 
from algebraic to exponential. We thus conjecture that 
while pinning generates metastable states, the internal 
organizations of the defects with respect to a mesoscop- 
ically ordered state give rise to an exponentially large 
number of these metastable states. An extensive S c does 
not occur in the case of pinned defects about a micro- 
scopic lattice because there are no internal degrees of 
freedom. 

In summary, we have demonstrated glassiness (in the 
sense of a finite configurational entropy) in a system in- 
teracting through a Kac-type repulsive step-potential. 
We show that although the mesoscopic modulation 
length Iq is due to the interaction potential alone, in con- 
trast to the model proposed in Ref. M, the underlying 
principle for glass formation is identical, i.e. £ > 21q. 
In the absence of a short-range interaction (the gradi- 
ent term in the Hamiltonian) this criterion for glassiness 
can only be fulfilled when V(q) possesses at least one 
minimum. The concept of a defect-wandering length is 
useful for understanding the source of the exponentially 
large number of metastable states. The picture of the 
origin of the configurational entropy as the organizations 



of defects in ordered arrays of mesoscopic structures is 
consistent with the numerical findings in Ref. fTij . The 
glassiness we have found is due to configurations of meso- 
scopic defects. Based on the analysis of a ID Ising model, 
it is speculated that only pinning of defects about the 
mesoscopic pattern can give rise to exponentially large 
numbers of metastable states. 
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